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Abstract 

We consider the Schrodinger operator on the zigzag and armchair nanotubes (tight- 
binding models) in a uniform magnetic field and in an external periodic electric 
potential. The magnetic and electric fields are parallel to the axis of the nanotube. We 
show that this operator is unitarily equivalent to the finite orthogonal sum of Jacobi 
operators. We describe all spectral bands and all eigenvalues (with infinite multiplicity, 
i.e., fiat bands). Moreover, we determine the asymptotics of the spectral bands both 
for small and large potentials. We describe the spectrum as a function of \^\. For 
example, if — > ^(f ~ 7^ + 7rs)tan2^,A: = 1,2, ..,A'", s G Z, then some spectral 
band for zigzag nanotube shrinks into a fiat band and the corresponding asymptotics 
are determined. 



1 Introduction. 

After their discovery [li], carbon nanotubes remain in both theoretical and applied research 
[SDDj . Structure of nanotubes are formed by rolling up a graphene sheet into a cylinder. 
Such nanomodels were introduced by Pauling |Pa] in 1936 to simulate aromatic molecules. 
They were described in more detail by Ruedenberg and Scherr [RSlj in 1953. Various physical 
properties of carbon nanotubes can be found in [SDDj . 

There are mathematical results about Schrodinger operators on carbon nanotubes (zigzag, 
armchair and chiral) (see |BKj . |KLj . [KLlj . |Klj . [KuPj . |Pkj ) . All these papers consider the 
so called continuous models. But in the physical literature the most commonly used model 
is the tight-binding model. ("In solid state physics, the tight binding model is an approach 
to the electronic band structure from the atomic limit case. In the tight binding model, it 
is assumed that the Fourier transform of the Bloch function can be approximated by the 
Linear Combination of Atomic Orbital(LCAO). Starting from the Hamiltonian of an isolated 
atom centered at each lattice point, the band structure of solids can be investigated.") For 
applications of our models see ref. in [ARZj . [SDDj . |Haj . 
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In this paper we concentrate on carbon nanotubes which arise from graphene: zigzag 
and armchair nanotubes (see physical propereties in (SDDj ). We will study and compare 



spectral properties of Shrodinger operator on zigzag and armchair nanotubes. We will show 
that these operators have different spectral properties. 
For example: 

1) The Shrodinger operator Hzi on the zigzag nanotube is unitarily equivalent to the 
direct sum of scalar Jacobi matrices (see Theorem 12. ip . But the Shrodinger operator on 
armchair nanotube Har is unitarily equivalent to the direct sum of Jacobi matrices with 
2x2 matrix valued coefficients (see Theorem 16. II) . Then the spectral analysis of Har is more 
difficult. 

2) For some amplitude of the constant magnetic field the spectrum of Hzi has absolutely 
continuous part and eigenvalues (flat bands, see Theorem 12.21) . But the spectrum of Har is 
purely absolutely continuous for any amplitude of the magnetic field. 

3) The spectral bands of operators H^i and Har are different. But in some cases the 
spectra of these operators has the same part (see Theorem 16. 2p . 

4) In the simple case, when the magnetic field is absent and external electric potential has 
minimal period 2 the spectrum of H^i and Har are coincide. Remark that the multiplicity 
of some spectral zones is different (see Sect 4 and Sect. 6.2). 

5) The structure of spectral zones of Har and Hzi for large electric potentials is similar, 
since the spectrum is a union of small clusters, but asymptotics of this clusters are different 
(see Theorem 12 . 6 1 and Theorem 16. 51) . Moreover, we have similar situation for small potentials. 

In the proof of our theorems we determine various asymptotics for periodic Jacobi op- 
erators with specific coefficients see f l2.5l) . Note that there exist a lot of papers devoted to 
asymptotics and estimates both for periodic Jacobi operators and Schrodinger operators see 
e.g. IKKulj . [Li], [OTu], [SI], [S2]. 



2 Zigzag nanotube. 

In this Section we consider the Schrodinger operator H^ on the zigzag nanotube F C M'^ (ID 
models tight-binding model of zigzag single- wall nanotubes, see |SDD] . [N]) in a uniform 
magnetic field = \^\eo, Bq = (0,0,1) G and in an external electric potential. Our 
model nanotube F is a graph (see Fig. [2] and 2) embedded in oriented in the 2;-direction 
eo with unit edge length. F is a set of vertices (atoms) connecting by bonds (edges) Tn^kj 
and 

r = U^e^r^^, rn,o,fc = Xn+2fc + -^eo, r^.i.fc = r^^cfc + eo, uj={nJ,k)eZ, 
Z = Zx {0,l}xZ^, Zn = Z/{NZ), xfc = i?(cos^,sin^,0), i? = (2.1) 

Our carbon model nanotube is the honeycomb lattice of a graphene sheet rolled into a 
cylinder. This nanotube F has hexagons around the cylinder embedded in M^. Here 
G Z labels the position in the axial direction of the tube, j = 0, 1 is a label for the two 
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Fig 1. A piece of zigzag nanotube. 



types of vertices (atoms) (see Fig. [2]), and k E labels the position around the cylinder. 
The points ro,i,fc, k G Z^v are vertices of the regular N-gon and ri o,fc are the vertices of 
the regular N-gon J^i. arises from by combination of the rotation around the axis 
of the cylinder C by the angle and of the translation by |eo. Repeating this procedure we 
obtain F. 

Introduce the Hilbert space £^(F) of functions / = {fuj)ujGZ on F equipped with the norm 
Il/llf2(r) ~ X]wg2 l/i^P- The tight-binding Hamiltonian on the nanotube F has the form 
H'' = Hq + V on ^^(F), where Hq is the Hamiltonian of the nanotube in the magnetic field 
and is given by 

(-f^o/)™,0,fe = e.^^^ fn-l,l,k + e*^^/n-l,l,fc-l + G^'^^ fn,l,k-i 

(-f^o/)n,l,fc — ^ *^Vn+l,0,fc+l + e *''Vn+l,0,A: + 6 ^^'^ fn,0,k, f = ifui)u)i^Z, 

uj = {n,j,k) eZx {0,1}xZn, h = 0, 61 = -62 = 6 = ^ cot (2.2) 

lb ziV 

(the last line in (12.21) was obtained in [KLl ]) and the operator V corresponding to the external 
electric potential is given by 

(Vf)^ = V^f^, where = ^2™, K,o,fc = '"sn+i, v = {vn)n& e (2.3) 

Such potentials can be realized using optical methods, by gating, or by an acoustic field (see 
[N]). For example, if an external potential is given by Aqcos{C,oZ + /?o) for some constant 
Aq, ^0, Po, then we obtain 

V2n = A cos(2<(n V2n+i = A cos(27r^n + p), n G Z, (2.4) 

for some constant A, If ^ is rational, then the sequence Vn,n E Z is periodic. If ^ is 
irrational, then the sequence Vn,n eZ is almost periodic. 
Below we use notation for the set {1, .., j}, j ^ 1. 
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Fig 2. Nanotube in the magnetic field. 



Theorem 2.1. i) Let v = {vn)n&'L ^ Then each operator , b is unitarily equivalent 
to the operator ®i J^, where is a Jacobi operator, acting on and given by 



{Jky)n = an^lVn-l + anVn+l + VnVn, U = {Vn) 



0-2n = 0-k,2n — 2 Cfc 



a2n+i = ak,2n+i = 1, Cfc = cos(6 + 



iV' 



neZ, (2.5) 



"^k+i' ~ ^N-k /^'^ '^^^ {k,b) G Z^v X R. Moreover, the operators and 



and jI^^ - 

H'^'^n are unitarily equivalent for all 6 G M. 

irk 
N 



a) Let, in addition, Ck = cos(6 + ^) = for some (k, b) G Z^r x M. Then 



<Jl) = ^PpiJk) =|^n,, = + i-iy\vf + 1|^ - 
Remark. 1) The matrix of the operator J| is given by 



V2n-1 ± V2r^ 



(2.6) 



J I. 
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(2.7) 



2) If |cfc| = |, then is the Schrodinger operator with a„ = 1 for all n G Z. In particular, 
if 6 = 0, ^ G N, then J% is the Schrodinger operator. 

3) In the continuous models similar results were obtained in |KLj . |KLlj . 

4) Exner [Exj obtained a duality between Schrodinger operators on graphs and certain Jacobi 
matrices, which depend on energy. In our case the Jacobi matrices do not depend on energy. 
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1. Periodic electric potentials v. Introduce the class P'J^'^ of real s-periodic sequences 
matrix where 



V = {vn)nez £ and Vn+s = Vn, for all n E Z. If v G t^^J-,p* ^ 1, then J| is 2p-periodic 



p={'^ ^/ ^3 even ^^.g) 
p* if p* is odd 

If Cfc 7^ for some {k,b) E Ijn >^ then the spectrum of J| has the form 

see |vM] . where z^'^ are 4p-periodic eigenvalues for the equation a„_iy„_i + a„?/„+i + f„y„ = 
2/ = (2/n)nGZ- The intervals are separated by a gap 7^,, = {z\~^,z\'^^) of 

length ^| ^0. If a gap 7^ „ is degenerate, i.e., [7^ „| = 0, then the corresponding segments 

If Cfc = for some (A;, h) G Zjv x M, then (12. 6p gives (t(J|) = cTpp(J^), where 

M4)={^n,,=< + (-iy\A?7I, = !^?!izi±^, (n,j)eN,xN2|, (2.10) 

and each eigenvalue of J\ is a flat band, i.e. it has infinite multiplicity. In Theorem 12.21 
we show that the spectral band shrinks to the fiat band {A„} as — > and the 
corresponding asymptotics are determined. 

Each operator J\ is unitarily equivalent to the operator I'^^-k) -^l^**) = 2|cfc|, where 

2p X 2p matrix K{t^ a) = K{t, a, v) is a Jacobi operator, acting on C^^ and given by 

/ 1 ... f \ 
1 a ... 
a ... 

\ ra ... 1 / 



i^(r,a) = ir°(r,a) + 5, ir°(r,a) 



5 = diag(t;„)^ (2.11 



where r G S"*^ = {r G C : |r| = 1}. Let /ii(r, a) ^ fi2{T,a) ^ ii^{T,a) ^ .... ^ /i2p(T, a) 
be eigenvalues of K{T,a),T G §\ here fin{-,ci) is analytic function in r G §^ Note that 
/i„(S^,a) = cr^n for all (fc,n) G Zat x N2p. If Ck 7^ 0, then each /i„(-,a),n G is not a 
constant and |(T^ ^| > 0. If = for some k G Zjv, then each /!„(■, 0) = const = A„, n G 
and cr^ „ = {A„} is a flat band. 

2. The case = 0. Consider the Schrodinger operator at ^ = 0. By Theorem 12.11 
the operator is unitarily equivalent to the operator ®i J^, where J° is a Jacobi operator 

at 6 = and here a2n = 2| cos^|, a2n+i = 1- Note that if 7^ y, then cr(J^) = aadJ^) 
and if k = then (t(J°) = o-pp{J^). 

3. Example of simple periodic potentials t;. Consider the potential v = V2k+i = — ^2fe ^ 
M, /c G Z. In Section 4 we will show that 
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4J = ±v/^2 + (2|c,| + l)2, = ±^v^ + {2\cu\-l)\ k G Z^, 

where 7^^^ is the gap in the spectrum of J|. This gives 

[|±Vl+^^^| Cfc = 0, 
and then we deduce that the spectrum of is given by 

am = a^AH-, U = (» ^' ^ °- ""^ , (2.12) 

I |±v 1 + "i^ I ^/ Ck = U, some fc G Z^r 

= 4n \ 7(if^), 7(^') = (4'-, 4^"-), (2.13) 

where '^{H^) is the gap in the spectrum of . Note that if = for some k G Z^v then 
= {±yrT^} C i{H^). Theorem Oi yields a{H^+^ = a{H'') for all 6 G M. 
Then we need to consider only the case 6 G [0, and in this case we get 



l-J ^ 2iV 



4^=rK L (2-14) 



Moreover, in particular case = we obtain 

-i[H^) = {-\v\,\v\). yGN, 6 = 0. (2.15) 

Now we return to the general case of periodic potentials. First theorem is devoted to the 
asymptotics of small spectral bands that degenerate to the flat band. 

Theorem 2.2. Let v G i^l"^ , ^ 1 and Ck ^ as b ^ Bq = ^ — ^ for some k G Z^r and 

let Ai ^ A2 ^ .. ^ be eigenvalues of K{l,0,v). Then the endpoints zl'~ , s G N2p 

of the spectral bands a^ ,, = zl'~] are analytic functions in b & {\b — bo\ < e} for some 

e > and satisfy 

%::_! = A. + 0(4), z',- = \s + Oicl) as Cfe^O. (2.16) 
Let in addition Xs-i < Xs < K+i for some s G N2p, where Aq = —00, A2p+i = +00. Then 
2 

4':: = - j^\2ck\^ + C^M(2cfe)'" + 0«+^), A,= n (2-17) 

* 2^2n^p n=l, n^s 



\<s\ = <; - 4U = + (2.18) 

as Ck ^ for some constants Ck^n, which depend only on v. 
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Remark. By fl2.16p . each spectral band crl^,n E shrinks to the flat band {A„} as 
Cfc 0. 

We consider the nanotube in weak electric fields. Our operator has the form H^{t) = 
Hq + tV , where a coupling constant t — 0. In this case the corresponding Jacobi operator 
depend on t and is given by 

{Jl{t)y)n = an-lVn-l + a-nVn+l + tVnVn, V = (yn)nGZ ^ £^ n G Z, (2.19) 

and a2n = 2|cfc|, a2n+i = 1- We study how the spectral bands cr^ „(t) = [^^ n-il^)' ^t' n (^)]? 
n e of the operator J|(t) depend on the couple constant t ^ 0. 

For V G i^l"^ we define two vectors f ° = (f2n)i, = ('i^2n-i)i ^ and 

M„ = (M,e„), M G C^, e„ = ^(r^^)^^;^ G C^, r„ = e*"?, Up+n = Up-n, n e Np. (2.20) 

Here e„, n G Np is a basis in and {u, e„) is the scalar product in C^. Define -^qX = {v E 

N2p_i tf 2\ck\ = 1 

N2p-i\M tf2\ck\^l' 

Theorem 2.3. Let Ck ^ for some {k, 6) G Zat x R. Let v G i^^^p^ and let f ° = (f2n)?, = 
('^2n-i)i £ H^^- r/ien i/ie asymptotic of the spectral hands a\^{t) = [%.'^_i(^), ^.'^(i)]? € ^2p 
of the operator J^{t) hold true 



'pT ■ Ef = 0} and the sets Nk,p 



&) = <.(0) ± tAA^) + 0{t'), n G Nfc, 



If.O _|_ „2iarg(2|cfc|+T„) -1 1 ri ^ n 



<l\, 2|cfc| = 1, n=p 



< W = < (0) + O(t^), (t) = 42p(0) + O(t^), 

and ^/2|c,|^l ^ 4JW = 4:^(0) +0(t2), (2.22) 

4:^(0) = |2|cfc|+r„|sign(n-p), n G N2p-i \ M, zj;j(0) = ±|2|cfc| - 1|, (2.23) 
as t [Q. Moreover, if p^ is odd, then for all n G fik,p the following identities hold true 

/|(-l)V„ + e--(^l-l+-)|, n^p 1^'=-^°' 

Pfc.n = < ..oi I 1 J ,SPk,ntO, if\ck\ = ^, even n . (2.24) 

10, tj 2\ck\ = I ana n = p i i i 

[pk,n = 0, if \ck\ = 2, odd n 

To describe some examples of external fields which create the open gaps we define the set 

X ^( ^^per .h'^ + vl 7^0, vy^ = 0, aline Rp^,, {;0^0,p,G2N U 
I °'^*'|{;0^0, allneKp-i, p, is odd J J 
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Proposition 2.4. i) The set Xp^ 7^ for any ^ 2. 

ii) If V G Xp,,p* G 2N, then 

= 4,n{^)±t^n + 0{t^), Cn=\Vn + Vn\>0 as t 10, all H & Nk,p. (2.26) 

in) If V & Xp^ is sufficiently small and is odd, then 

If2\ck\ ^ 1; then each i)k,n{v) 7^ 0, n G N2p_i \ {p} and „ 7^ 0. 

If 2\ck\ = 1, then each i^kniv) = \ ^ ^ and '~fkn ^ for any even 

10 all odd n G N2p_i 

n G Nap-i. 

Remark. (12.261) gives the asymptotics of the gap length -2fc'^(t) — ^(^) = ^2|{)^| + O(t^) as 
t — s> where j = or j = 1. Note that the first term does not depend on A; G Z^v- If is 
even, then for large class of potentials f G all gaps (2;^'~(t), ^^^(t)) are open. 
We formulate the theorem, motivated by the physical paper of Novikov [N] . 

Theorem 2.5. Let v G l^^J and let t > 0,b be sufficiently small. 

i) Let 6 = 0. If N E 2N and p are coprime numbers, then app{H^{t)) C r\^^lcr{J^{t)). 

ii) If p > 2N , then the spectrum of H''{t) on the set a{H^{t)) fl ([— p, — r] U [r,p]) has 
multiplicity 2 and satisfies 

2p—l I 12L I 

a{H\t))n[r,p]=criJ'^{t))n[r,p] = [r,p]\ [j 7^^), r = |2 + |, p = ^±^±^, 



JL 

N 



a{H\t)) n [-P, -r] = a(4(t)) n [-p, -r] = [-p, -r] \ (Jt^.J^)- (2.27) 

1 

Moreover, if v E Xp,, t/ien each |7Ar,n(^)| > 0,?7, G N2p_i. 
Hi) IfN^ m, then cr{H^{t)) n [-r,r] = for some r > 0. 

iv) If N E 3N anc? p > 2N, then the spectrum of H^{t) on the set a{H^{t)) fl [— r, r] has 
multiplicity 2 and satisfies 

a{H\t))n[-r,r]=a{Ji{t))r\[-r,r] = [-T,r]\ (J 7|,„(t), r=|l-et|, (2.28) 
I r+M^n / P*e2N, nGN2p-i 

l74,nWI>0 ^/ S . > veXp,. (2.29) 
3 I odd, even n G N2P-1 

Remark. 1) The gaps 7Ar„(t) in (12.271) and 7^ {t) in (I2.28P are also the gaps in the 

' 3 

spectrum of H^{t). Then we may choose the potentials v such that all these gaps are open 
(for wide set of potentials). 2) Due to iii) cr{H^) has a gap contained the interval [— r, r] 

We consider the nanotube in strong electric fields. Our operator has the form H'^lt) = 
HQ + tV, where a coupling constant t — 00. For each (fn)i^ ^ 1^^^ there exists a permutation 



8 



a : — >■ such that /i„ = and hi ^ h2 ^ ■■■ ^ h2p. Let f„ 7^ for all n 7^ j, 
n,j G N2p. Defining disjoint intervals ^„ = /;,° = ^"+^"+1 , n e Nap, = 

— 00, /;,2p+i = 00, we obtain the inclusion a{H'^{t)) C U^^q^„ = R. We shall call the set 
a{H^{t))r\^n the n'th spectral hands cluster. Our goal is to study the asymptotic distribution 
of eigenvalues in the n'th cluster as t ^ 00. 

Theorem 2.6. Let v E i^^J' , f„ 7^ vj for all n 7^ j , n,j E and let Ck = cos(6+ ^) 7^ for 
some {k, b) E N^f x M. Let Va{n) < '^a{j) for all n < j and some permutation a : — > Nap. 
Ifn = a^^{n) for some n E Nap, then the spectral hands cr^„(t) = [-2fcn_i(^)) z'k^{t)] satisfy 

4:l-.W=t''~- ^" c- = ;7^Hr + ;7^. (2-30) 

w - -K-iW = E- = ^,\{ - (2-31) 



as t ^ 00. Moreover, 

N 



a{H\t))n^n{t) = U^'nW C U„t--,t^„t+- , 5 = max- -, (2.32) 

[|cfc| 7^ Icfc'l, oE^N 

where the spectrum of H^{t) on cr\~{t) has multiplicity 2 if Ck and (y\~{t) is a flat hand 
ifck = 0. 

Remark. 1) Theorems 12. 3^ 1.4 describe the case t and Theorem 12.61 describe the 
case t 00. These two cases are quite different, see Fig. 3 and Fig 4. 



-3 3 

Fig. 3. Open spectral small gaps for the potential tV as t ^ 0. 



Fig. 4. Spectral clasters for the potential tV as t ^ 00 for the case = 4. 



2) The spectral bands cluster a{J^(t)) fl ^n(t) is a union of non overlapping bands 
(T^~(t),A; E Njv, see fl232|) . Recall that if \ck\ = \ck'\, then Jl{t) = J^,{t). 

We present the plan of our paper. In Sect. 2 we prove Theorem 12.11 and 12. 2[ In the 
proof Theorem 12.11 we use arguments from [KLj . [KLl] . In the proof Theorem 12.21 we use 
arguments from |KKul] . In Sect. 3 we consider the simple examples for the case p = 1, in 
fact, we study unperturbed Hamiltonians. In Sect. 4 we prove Theorem 12.31 - I2l6l In Sect. 
6 we apply some of these methods to analyze the spectral properties of Shrodinger operator 
on armchair nanotubes. 
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3 Proof of Theorems 12.11 and 12.2 



Proof of Theorem l2.1l i) Define an operator : (f^)^ acting on a vector-valued 

function ^ = (V'n)nez e (^^)^,^2n+i = (/n,o,fc)fceZAr, V'2n = (/n-i,i,fc)fceZiv ^ by 



(3.1) 



Define a matrix-valued operators P„ : by 

P2n+lh = iVn,0,khk)keZMJ P2nh = (Ki,l,A:^fe) fcGZ]v ) h = {hk)k&N- (3-2) 



Define the operator S in C by Su = {u]y,ui, . . . ,un_i) , u = (Mn)i G C . Using (13. II) . 
(12:21) . (1I3D,(D and 5* = A = e'^I^ + e-'^S* we obtain 

(^V)2n+l = (e*^5 + e-^^)^/'2n + ^2n+2 + i^2n+lV'2n+l = ^>2n + V'2„+2 + i^2„+l V'2n+1 , 
(^V)2„ = V'2„-l + (e* + e-*5*)^2n+l + ^'2nV^2n = V'2n-1 + ^^2n+l + AnV^2n. 

Finally we rewrite the operator : ^ in the form of the operator Jacobi by 

(^V)n = A:-l^n-l + An^n+l + Pn^n, = A = + C-^'^* , ^2^+1 = Jtv, (3.3) 

n G Z, and then 



V 



A* Pi /;v 

Q In P2 A 

A* P3 In 

Jiv P4 

A* 



\ 



(3.4) 



The matrix- valued function P„ is 2p-periodic. Then the operator is a 2p-periodic Jacobi 
operator with N x N matrix -valued coefficients. Note that such operators were considered 
in |KKu2j . 

The unitary operator S has the form S = Yli ^^'^k-, where Se^ = s^'^k and = 
-^{1, , s~'^^ , s~^^^^Y is an eigenvector (recall that s = e^^)] VkU = ek{u,ek),u = 

{un)i G is a projector. Define the operators Sip = {S'ipn)n& and Vki' = {'Pk'ipn)n&- 
The operators S and commute, then = ®\{^^Vk)- Using (13. 3p . (13.41) we deduce 
that ^^Vk is unitarily equivalent to the operator given by 

Ck = cos(6 + ^), ak,2n = e* + e'^^s-^ = 2e-^^Cfc, s = e'^ , ak,2n+i = 1, (3.5) 



and using Lemma [3TT] we obtain (12.51) . 
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ii) If Ck 



0, then the Jacobi operator has the form 






Vl 1 











1 V2 
















1 








1 

















V 



The eigenvalues of Jn are given by Znj 
Z X N2, which yields ^Mj- ■ 



V2n-1 1 
1 V2r, 



(3.6) 



Recall results from [vMj about our 2p-periodic Jacobi operator J(c 



for (n,j) G 
^ given by 



{J {.a)y)n = an-lVn-l + anVn+l + VnVn, = CL > 0, Og. 



n+l 



1, neZ, y = {yn)n&- (3-7) 



Note that = J{a), where a = 2\ck\,Ck = cos(^ + b). Introduce fundamental solutions 

(f = {!fn{z, a))n& and d = {^n{z, a))nez for the equation 



an-lVn-l + anUn+l + VnUn 



{z,n) G C X Z, a2n+l = 1, a2n 



(3.8) 



with initial conditions ipo = 'di = 0, ipi = {}o = 1. The function A = ^{ip2p+i + ^2p) is 
called the Lyapunov function for the operator J(a). The functions A, y9„ and 'i^n,''^ ^ 1 are 
polynomials of {z,a,v) G C'^^'^'^. It is well known that a{J{a)) = (Tac{J{a)), where 



cXaciJia)) = {z eR: A{z,a) G [-1,1]} = 



UiV„(a), 



[Zn~M),z^{a)], (3.9) 



and Zq < Zi ^ z^ < .. ^ Z2p, where z^ = z^{a). Note that A{z^,a) = (— 1)p~" for all 
n = 0, ..,p. Below we will sometimes write a{a,v), J{a,v), .., instead of a(a), J(a), .., when 
several potentials are being dealt with. Recall that the 2p x 2p matrix K{t, a) is given by 



K{T,a) = K\T,a) + B, K^r^a) 



/ 1 ... f \ 

1 a ... 

a ... 

\ ra ... 1 / 



B = diag(^;0?^ (3.10) 



where r G = {r G C : |r| = 1}. Fix a,(j) ^ [0,27r], then eigenvalues of K{e'^'^,a) are 
all zeros of the polynomial A{z,a) — cos (p. Then the fundamental solutions (pk,n,'&k,n, the 
Lyapunov function and the spectral bands for the operator J| satisfy (see also fl2.9p ) 



'fk,n = 'fn{z, a), 'dk,n = i^niz, a), Afc = A{z, a) z^;^ = z^{a), (3.11) 
<^iJt) = cTaciJt) = {zeR: A,{z) G [-1, 1]} = at,,, a^,, = [^^i, 4';], (3-12) 
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Proof of Theorem 12.21 Let a = 2\ck\ — »• 0. We consider the matrix K{t, a) as a — > 0, r G 
§1 = {r G C : |r| = 1}. If a = 0, then we get K{t, 0) = ©?X, where X is given by ([XSD- Let 
-^1 ^ ^2 ^ ••••A2p be the eigenvalues of K{t, 0). The endpoints z'^_^{a), z~{a) of the spectral 
bands (T„(a) = [z^_i{a), z~{a)] of the operator J(a) are the eigenvalues of i^(±l,a,t>). By 
the perturbation theory |RSj . they are analytic function from a and if a ^ 0, then the 
spectral bands converge to the set {Ai, A2, ....A2p}. The number of spectral bands converging 
to \n coincides with the multiplicity of A„ as a — 0. In particular, if some A„,n G is 
simple, then (T„(a) {A„}. 

Recall that the monodromy matrix for the operator J(a) is given by 

M.A^) = (/^^ ) = T,.T2T„ 

Tn=-{\ " V° M = f ~" V (3-13) 

a\-l Z-V2n+lJ \-a Z-V2nJ \ V2n - Z (pn/ a J 

where 0„ = - V2n){z - f2n-i) - 1- Let 



Then M2p = E~^XpAXp_iA..AXiE^^ , which yields the Lyapunov function A given by 

1 1 ^ 

2A = TrM2p = TTXpAXp_,A..AX,A = -Ti XpA,Xp-iA,..AX,A, = _^a2"$„(^) 



ra=0 

and 

2a^ 



A(z,a) = ^°^"^ + "''^^'"\ $o = n(^-^'^)' $(^,t) = 5:r-$.(.), (3.14) 



n=l n=l 



for some polynomials $„. By the perturbation theory (see |RS] ) . the endpoints z^,Z- of 
the spectral band crs(a) = [zf_i{a) , z~ (a)] = are analytic functions in some disk 

{a G C : |a| < e},e > and satisfy the equation A{z±,a) = =]=(— 1)*, which has the form 

%{z±) + a2$(z±, a^) = T{-iy2aP. (3.15) 

Moreover, they satisfy z±{a) = A^ + O(a^) as a — ^ at p ^ 2 (see the case p = 1 in Sect. 3). 
Let As be a simple eigenvalue for some s G N2p. The differentiation of ( 13.15^ yields 

z'^{a)n + da{a^^{z±, a^)) = Tp(-l)'2aP^\ n{z, a) = d,{%{z) + a^^{z, a^)). (3.16) 

The differentiation of fl3.15p r G [l,p] times yields 

z^^^n{z^) + Gr{zt'\ a) = ^^-ZL_(_i)^2a^-, (3.I7) 
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for some polynomial Gr- Then at a = this gives 



{r-l), 



pl 



(p — r)! 



-iy2aP- 



\a=0- 



(3.18) 



Thus we obtain ^^^'^^^''(0) = for all 2r+l < p, since the polynomial $ = ^{z, a^). Moreover, 
using Z-{0) = -2+(0) we obtain z^l\o) = z^\o) for all r < p. 
Consider the case r = p. Identity fl3.17p implies 



z^^\0){-irAs + G,{zr'>{0),...,z+{0),0) = Tpl2{-iy 



(3.19) 



which yields z^\o) = j-{Cp =F 2) for some constant Gp G M. Using this and as{a) = 
[z+_i(a), z-{a)] = [z+, zJ\ and fl3TTD we obtain ^TH) . ^TM . ■ 

Lemma 3.1. Let a Jacobi operator J : l"^ ^ (? is given by 

{Jy)n = a*„i?/„_i + anVn+l + VnUn, V = {yn)n& & ^^ = a„ G C, f„ G M, (3.20) 

n G Z, for some p ^ 1. Then 

= J^, iJ^y)n = \an-l\yn-l + knll/n+l + VnVn, (3.21) 

where the unitary diagonal operator \1' is given by 

% = (Mnl/n)neZ, = _[ | , « ^ 0, Un = [[ej , U < , En = V^"^ „ • (3-22) 



if an = 



Proof. Direct calculations give (13.211) . 



4 Example for the case p = 1 

In this section we consider the Jacobi operator J^^k ^Ijn given by 



t-'i. 



v 



a 1 

Q I -V a 

a 1 

1 -f 

a 



a = 2|cfc|, w = f2n+i = -W2n e M, ?2 G Z, (4.1) 



i.e., the case p = 1. The monodromy matrix M2 satisfies (see f l3.13p ) 
f2\ f 1 \ f 1 



M2{z) 



1^3 



1 

1_ z+v 

a a 



—a Z — V 



—a 
—z — V 



z — V 



(4.2) 
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Let A*^ = = ^ be the Lyapunov function for the case a = 1. This yields 



TrM2 + 4 ,^ ttA:, ,^ ^, 

^^--(^+iv)- ^^-^^ 

The periodic eigenvalues z^'q satisfy the equation Ak{z) = 1 and anti-periodic eigenvalues 
z^'^ satisfy the equation Ak{z) = —1 and they are given by 

zl'^ = ±^v^ + {2\c,\ + ir, z]:^ = ±./^^^TWJ^. (4.4) 
The spectrum of J| has the form 

'^(Jk) = [4fl^ u [zl'^, 4;o] = [4;o ' \ 7^,1, (4.5) 

where 7^ ^^ = (-Zfc'i '-^fct) ^ S^P- Note that 

7ti = (4:r,^fc:t)^0, ^/ (4.6) 

Let Cfc ^ 0. Then (12T7D . ( IZTSi) yield 

4|cfe| 



\4 \ = - = — — + 0(4), w = vi 



w 



^fe 1 = 



+ ^ + 0(4), .tJ = -^-^ + 0(4). (4.7) 

w ' w 

1. The operator H^, no magnetic field, 6 = 0. In this case using (14. 4p . (14. 5p . we obtain 

^0 z/A:^{f,f} 



0,+ 0,+ 0,± _ 1/2 , q ^, , _ f/.- /.+^ J - ^ ^ 1 3 ' 3 J (Ai 



and then 



{(/) f :^ (3i 1^ 

a..(i7°) = [.oV,<o1\7(i^°), 7(^°) = l%,_ 0,+ ,^, Zl^^^ (4-10) 



N 



for some m G Z^v, where roughly speaking m - ^ 
2. Magnetic field, 6 7^ 0. Using i) of Theorem O we obtain a{H''+^) = a{H^), beR. 
Then we need to consider only the case b G (0, ■^). Using fLT\ we obtain 

.m = a^iH^)<Ja„iH% .„iH^) J' t''„'""' fj, .(4-11) 

I |±v 1 + V'^^ I J Ck T some k G Zat 

14 



a^^{H') = [z',\ zj'-] \ ^{H'), j{H') = {z\'~, z?'"^), (4.12) 
where 'y{H^) is the gap in the spectrum of and 

1j , (4.13) 
and 

-i{H') = {z'{-,z'{-'), z\^ = ±V^2 + (2|cfc| - 1)2, /or some A; G Z^, (4.14) 
where roughly speaking 2\ck\ ~ 1. 



5 Proof of Theorems [SSHSIEl 



Proof of Theorem 12.31 In order to determine the asymptotics fl2.2ip we need the following 
fact from the perturbation theory [RS] : Let A{t) = Ao + tAi, t G M, where Aq = Aq, Ai = Al 
are operators in C^^. Let /i be an eigenvalue of Aq of multiplicity 2 and let be the 
corresponding orthonormalized eigenvectors. Then there are 2 functions /i±(t) analytic in a 
neighborhood of 0, which are all the eigenvalues. Moreover, /i±(t) = fi + /i±(0)t + O^t"^) as 
t — >■ 0, where n'±{0) are the eigenvalues of P*AiP and P = {h~ , h~^) is the 2p x 2 matrix. 

We determine the asymptotics (12.211) of zl'^{t) for k e Np,n ^ 0,p,2p, the proof of 
other cases is similar. We apply the perturbation theory to the operator K{±l,a,tv) = 
K^{T,a) + tB as t ^ 0, where K is given by fl3.10p and a = 2\ck\- Recall that z^^{t) 
are eigenvalues of K{±.l,a,tv), (see (13. 7p - (13. lip ). The operators K^{±l,a) has eigenvalues 
^fcn(O) ~ ^fcn(O) ~ ^ti.'^) (with multiplicity 2) and the corresponding eigenvectors 

Zt = Znia), n e Z2,_i, (5.1) 

see Corollary 17.21 and (13. lip . Then by this fact, the derivatives (z^nYi^) eigenvalues of 
the 2 X 2-matrix P^^^BPk^n, where Pk^n = {Z'^^, Z^^) is the p x 2-matrix. Define the vectors 

F„ = (2p)-^(/,)f , /2,+i = r2V^-^(^l^'=l+-\ f2j=r'^\ r„ = eT, jeN,. (5.2) 

Let Vn = {v,Fn), n G Np and Vp+n = Vp-m n G Np_i. Using (15.11) . Corollary 17.21 (17. 6p we 
obtain 





^ 


Vn \ 




V Vn 


) 



Pk,nBPk,n=( ^^/"p^ ) = ( ^ )' ^^^^^ B = dmgivjf/ . 



TrB {b,Fn) 
{Fn,b) TiB 

The eigenvalues of the last matrix have the form ±|t;„|, which yields {z^j^)'{0) = ±\vn\- 
Recall that the orthogonal basis in is given by e„ = ■^{Tn)%i^ ^ where r„ = e~ 
and the vectors = (f2n)n=i and = {v2n~i)n=i^ ^n — (v-'^^n), n G Np, j = 0,1. Then 
flOl gives Vn = v^ + e2^^'-g('^+^"){;i and we obtain (K21^ . 
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Let S{ui, ..,Up) = (wp, Ml, .., Mp_i) be a shift operator. If is odd, then p = p^, and 
= S~v^ and using fl2.20p . we obtain 



since = r„ ^. Then we get 



2\ck\ = 1, then 1 + rP-ie2*^'-s(2|cfe|+r„) 



Simple calculations gives: if 2\ck\ ^ 1 and n G Nk,p, then 1 + T-P"^e^*'^''^(^''"=l+^"^ ^ 0, and if 

^ n is even 
0, is odd 

Proof of Proposition 12.41 i) Consider the case is even. Denote z = (zn)i for z = 
{znYi e CP. Using dOHD, we obtain ep_„ = e;;, n G Np_i and Cp = (2j9)-i(l, .., 1)^ g M^. 
If -u^ = X]n=i"«^« + "p^P' ^ = '^p-n ^ 0, n e Np_i, 7^ ftp G M, then G 
and -0^ = a„ 7^ 0, n G Np, since {e„}i is orthogonal basis in C^. Consider v'^ = —Vpep, 
then t>° G MP, since Cp G M'' and Wp = Op G M. Also -Up = — ttp 7^ 0. Then the vector 
V = vl, .., w^) G Xp,, since En=i(^n + ^n) = + = 0. Then Xp, 7^ 0. The proof of 
the case of odd is similar. The statements ii) and iii) follows from Theorem 12.31 fl2.24p . ■ 
Proof of Theorem EUl i) Using (El), we obtain that {1 - 5,1 + 5) U {-1 - 5, -1 + 5) C 
a(J^(0)) for any k G Nn-i \ {f } and for some 6 > 0. If A; = f then we obtain a{J^{t)) = 
app{H\t). Moreover, we have that app{H\t)) e {{1 - 6,1 + 5) U (-1 - S, -1 + 5)) for small 
t, then in order to prove i) we have to show that there are no gaps in small neighborhood of 
{±1}, i.e. we need to show that 2;^'^ {il}, 



|2|cfc| ^ 1, neNj 



p 



or 



or 



fevr nvr nvr 

2 cos h cos h I sm — 

N p p 



kn ( mr kir . 

4 cos — cos h cos — 7^ 1 

N \ p N ' 



rm k-K 

cos h cos— 7^0, (5.3) 

p N 

since cos 7^ for k ^ ^. The identity (15. 3p holds true, since p and are coprime. 

ii) Consider the case a{H^{t)) fl [— p, — r] the proof of other cases is similar. Theorem 12.31 
gives 

<Jl{^)) = [-2|cfc| - 1, -|2|c,.| - 1|] U [|2|cfc| - l|,2|cfc| + 1], ke Zn, (5.4) 
which yields (t(J^(0)) = [-3, -1] U [1,3] and 



[-p-5,-r + 5] C J^(0), [-P - 5, -r + 5] n J°(0) = 0, keN 



N-l 



for some small 5 > (see (12.51) and before (12.271) ). Then the spectrum in a{J^{t)) fl [— p, — r] 
has multiplicity 2 for all sufficiently small t and b. Also, using (12.231) . we obtain 2^„(0) G 
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Xn{t) =t{vn + eun.n + <yne^ + 0{e^)), = ""'^^^'" ^ lij^^, = (e°, K° (r)e° ^ , 



[— p, — r], 1 ^ ^ and z'^^i^O) ^ [—p,—r] for n > ^, which yield (12. 2 7p . The inequahty 
|7Ar,n(^)| > follows from Proposition 12.41 

iii) follows from (15 ■41) . since cr( J°(0)) fl [— r, r] = for any k and sufficiently small r > 0. 

The proof of iv) is similar to the proof of ii). ■ 
Proof of Theorem [2l6l Recall that {Jl{t)y)n = an- iVn-i + anVn+i + tVnVn, y = {yn)n& e 

n G Z, where a2„ = 2|cfc|, a2n+i = 1- Using fl3.10p we obtain 

Kk{T,a,tv) = t{B + eK'^{T,a)), B = diag{v as e = - — 0, a = 2\ck\. 
Then the perturbation theory |RS] for B + eK^{T, a) gives 

'Jj - Vn ' 

where i?e° = fje° and the vector e° = {5j^n)'^^=i G C^^. The definition of Uj^n yields 

Un-l,n = = C^n-l; l^n+l.n = 'Wn,n+1 = C^n; C^^*^ = if \j — n\ 1. 

These imply f l2.30p . since zl'^it) are eigenvalues of Kk{±l,a,tv). 

We show f l2.3ip for the case fi < .. < V2p, the proof of other cases is similar. Using the 
identity 2Ak{z,t) = 2A{z,a,tv) = det{zl2p — K{i,a,tv)) (see reasoning between (13. 7p 
and (I3.12P ). where w = YlT = |2cfc|%^^, we obtain 

2Ak{z, t) = det{zehp -B + eK\i, a)) = ^'^^^^'^^^^jp'^\ ^ = 7 = ^^'^^ 

where Fo(A) = det(A/2p-5) = ]\%i{\-Vj) and F is some polynomial of two variables A, e. 

Let A+(£:) = 2^'^_i(t)/t, A_(£:) = z'l^~{t)/t for some n G N2p. These \±{e) are the solutions 
of the equation F(A±, e) = ±1, where F(A, e) = Ak{z, t) and (l230D yields \±{e) = Vn + 0{e'^) 
as e — s> 0. By the perturbation theory |RSj, the functions \±{e) are analytic in some disk 
{|£:| < r},r > 0. Now we repeat the arguments from the proof of Theorem 12.21 after (13.151) . 
Differentiating (15. 5p 2p times we obtain 

(A+)(^)(0) = (A_)(^)(0), J < 2p, {Kf^\0) - (A„)(2^)(0) - 



En 



I.e., 

\X4s)-X4s)\ = — + 0is'^^') as s^O, E„ = J](^„ - ^;^.), 



p 



which yields (12.311) . since \+{e) = zl''^_-^(t)/t, \^{e) = zl'^{t)/t. Using (I2.30p . we obtain 
If |cfc| 7^ \ck'\, then (12.311) implies 

<„w n 4,„(t) = [z'^;n-iit),zi'-it)] n [4:'n-i(^),4:nW] = 

for sufficiently large t. This yields (12.331) for the second case. If k k' for k,k' G Njv and 
b ^ 2^N, then \ck\ ^ \ck'\ and we obtain (I2.33P for the first case. 

Using ([233D, we obtain a(4(t)) n (r{J^,{t)) = 0, k'. Then a{Jl{t)) has multiplicity 
2 and o"^_„(t) has multiplicity 2. ■ 

17 



6 Armchair nanotube. 




Fig 3. A piece of armchair nanotube. 




Fig 4. 3D model of armchair nanotube. 

We consider the Schrodinger operator H^{v) with a real periodic potential v on the 
armchair nanotube F C in a uniform magnetic field ^ = -8(0, 0, 1) G M^, 5 G M. Our 
model armchair nanotube F is a graph (see Fig. [6]) embedded in oriented in the ^-direction 
eg. F is a set of vertices (atoms) r^^ connecting by bonds (edges) and 

r = U^e2r^, cu = (n,j, A;) G Z = Zx {0,1} X Ziv, = Z/(iVZ), (6.1) 

where is a number of vertices in any ring of nanotube. The detail information about 3D 
coordinates of and about constants hj see in Appendix. 
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Introduce the Hilbert space i'^{T) of functions / = {fu})u}£Z on F equipped with the 
norm H/H^j^p^ = ^^^^ l/i^P- The tight-binding Hamiltonian H'^ (where b = (6i,62,&3)) on 
the nanotube F has the form H'^ = Hq + V on ^^(F), where Hq is the Hamiltonian of the 
nanotube in the magnetic field and is given by 

iH^f)n,l,k = e *''Vn+l,0,fc+l + e ^^'fn-lfi,k + 6 fn,0,k, f = (/a;)a;e2> 

u = {n,j,k) eZx {0,1} xZn (6.2) 
and the operator V corresponding to the external electric potential is given by 

(yf)^ = V^f^, where Vn,0,k = V2n, = W2n+1, keZN, V = {Vn)n& & 

(6.3) 

1. The operator H'' is an orthogonal sum of Jacobi operators. 

Theorem 6.1. Let v = {vn)n&z ^ Then the operator is unitarily equivalent to the 
operator ®i J\, where J\ is a Jacobi operator, acting on £^(Z) © ^^(Z) and given by 



02 ^ fJ-i 



{Jly)n = aVn-l + 0*yn+l + dnyn, y = {yn)neZ ^ 

a = ak= { Q l,s = eiv, rf„=( ^^^^^ 1 , n e Z. (6.4) 

Each Jji has absolutely continuous spectrum. 

Proof of Theorem 16. IL We give compressed Proof because this one is similar to the Proof 
of Theorem 12. 1[ Define the operator : (£^)^^ — ^ (£^)^^ acting on a vector- valued function 

^ = (^n)neZ e (iy, = {fnfl,k, fn,l,k)L^^ G C^^, by 

(^V)n = iiH'fU,k, (i^V)n,i,fc)Iez,. (6.5) 

Define the operator S in by Su = {u]\f,Ui, . . . ,M7v-i)^, u = {un)i G C^. Using (16.51) . 
(16.31) and S* = S^^ we obtain 

(^V)n = ^^n-l + v4*V'n+l + C„Vn, whcrC (6.6) 

. _ / e*''^^ \ „ _ / V2nlN e^'^N \ 7\ 



The unitary operator S has the form S = s'^Vk, where 

is an eigenvector (recall that s = e*i^); VkU = ek{u,ek),u = {un)i € is a projector. 
Define the operators Sip = {<Sipn)n& and Vk'ip = {^k'ipn)n&- The operators S and 
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commute, then = ®i {^^Vk)- Using (16.61) . (16.71) we deduce that ^^Vk is unitarily 
equivalent to the operator J^. ■ 

Below we use notation a = a{b,v) and dn = dn{b,v). 
2. The spectrum of unperturbed operator H^. 

We consider the case when all V2n+i = —V2n = v and 6 = 0, i.e. all Jk are 1-periodic 
Jacobi matrices. For this case we denote a = a{0,v), d = dn{0,v). The monodromy matrix 
for Jk is 



Mk(z) 



h ^ 

— (a)^ a{z — d) 



1 








1 





\ 













1 






_gk 





_gk 


{z — v)s^ 











z + v 


-1 


/ 



The determinant is 

Dk{z, t) = det{Mk{z) - Th) = / + t\s^ + 1) + Th\?> + ^ - z") + ts\s^ + 1) + s^"" 

= s-"^ (^^ + T^2Ck + -Z^)+ T2Ck + 1 

= s"^^ (^T + T-^f + 2cfc(r + r-^) + l + ^-z'' 

= s^^r^(r + r"^ - (2;))(r + r"^ - A^(2;)), where t = s~^t and 

A^(2;) = ±^ z^ -v^ - si- Cfc, where Ck = cos^, Sk = sm^. (6.8) 
The spectrum of Jk is 

a( J°) = G M : Dkiz,r) = for some r G S^} = (6.9) 
= {zeR: -2 ^ A^{z) ^ 2} = i-al) U (-a^ U (a^) U (a^), ^/lere 



^fc = ^5+^^+4^], al = [^v^ + sl,V5 + v^~Ack]. 

The spectrum of if is 



AT 



a{H) = [j a{Jk) = [-V9+^, V9+^] \ i-\v\, \v\). (6.10) 



k=l 

In particular case, ifv = 0, then <y{H^) = [—3, 3]. 

3. Small 2p-periodic real potentials. We consider the case 6 = 0. Firstly let J = J{q) 
£^(Z) is a p-periodic Shrodinger operator, i.e 

{JJ)n = fn-l + fn+l + Qufn, f = Un)n& ^ ^^(Z), 
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where q = {qn)n=i ^ ^rC^) 1n+p = In for all n G Z. It is well known (see |KKulj ). that 
the spectrum of this operator is absolutely continuous and has a form 

0"(J) = aac(J) = U?CT„, an = [Zn-l, z~], n G Np, (6.11) 

Zq < z]^ ^ < Z2 ^ < ... < z~. (6-12) 
We denote z^{q) = z^. Also we introduce spectral gaps 7„ = ■jn{q) as 

7„ = (2:^,2;+), neNp_i. (6.13) 

If g = then 

z^(0) = -2cos— , nENp-i, -4(0) = 2p-(0) =2. (6.14) 

For sufficiently small q we have (see |KKul] ) 

z±(g) = -2cos— + go± |gn| +0(||gf), g ^ 0, n e (6.15) 
P 



z+{q) = -2 + qo + 0{\\qf), z;{q) = 2 + qo + 0{\\q\\'), g ^ 0, (6.16) 

-2j\p-l 
n /j=0' 'n 



where we denote g„ = (g, e„), e„ = p H^n^)Ld5 i"n = e p 



Introduce the set Sp C M'' by 



{^a„(en + ep_n), all an ^0}. (6.17) 



Now we compare the spectrum of H^i{v) (zigzag) and H^j.{v) (armchair). 

Theorem 6.2. i) Let V2n = V2n+i, Vn+2p = for all n E Z. Let v^"" = (f2n)i o,nd J = J{v^^), 
then 

(a(J) + 1) U (a(J) - 1) C a{Hl). (6.18) 
a) Let N E 3Z and Vn+p = Vn for all n E Z. Let v = {vn)i and J = J{y), then 

a{J) C a{Hl). (6.19) 

Proof of Theorem 16.21 i) In our case (see (16. 4p ) we have 

4(0, «) = ( ' )=cf "V' \Ac- 

for some unitary matrix C {CC* = I2). Then (see (16. 4p ) unitarily equivalent to {J{v'^'") — 
I) © (J(f'^^) + /), where I is identity operator on P{Z). The statement ii) was proved in 
Theorem 12.31 (see also Remark 2) on page 5). ■ 

For example we describe the spectrum of H^{v) (armchair) near z = and near z = ±3 
for small potentials v (recall that a{H^{0)) = [—3,3]). 
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Theorem 6.3. Let V2n = V2n+i, Vn+2p = Vn for all n E 1^ and denote v^"" = (f2n)i- Let 
p > 2N > 4 and r± = 2 cos^| ~^ Jn ~^ ~ ^' ^hen for sufficiently small v we have 

a{H')n[r_,r^]= ([r.,r^]\ \J (7„ + 1)) U [[r_, r+] \ (J (7n-l)), (6.20) 

where first set and second set in the union has multiplicity 2. Also let 

- = l + 2cos(^ + l), f_ = l + 2cosl 

Then for sufficiently small v we have 

a(i7°)n[-f+,-F_] = [-f+,-F_]\ U (7n-l), (6.21) 



cr(/7°)n[F_,f+] = [f_,F+]\ U (7n + l), (6.22) 

where set on the right side has multiplicity 2. Moreover if v^'" G Hp then all |7„| ^ in 
(ICTl) - ( 16:221) . 

Proof of Theorem 16.31 We consider only the statement fl6.2ip . the proof of other 
statements is similar. We have (see fl6.14p ) 

-3 > -f+ > 2f(0) - 1 > -2^^1(0) - 1 > -F- > ^^jLi-^i(O) - 1 > (5 + 4ci)^ > -1. 

This inequalities shows (see (16.181) ) that for sufficiently small v we have 

a(i7°) n = a{J%) n -r_] = {a{J) - 1) n [-F^, (6.23) 

since n = 0, A; G N^v-i (see (El, O). Using identities ([S2SD and (ICT]) - 

fl6J[5|) we obtain fleA . ■ 

Let f be sufficiently small. We denote by Gar-, Gzi is a maximal possible number of 
the open gaps on the edge of spectrum, i.e. in the set a{H^^) fl [—3, —3 + a] and cr{H^^) fl 
[—3, —3 + a] respectively, where a is a some sufficiently small value. Now we estimate Gar-, 
Gzi for sufficiently large period 2p. 

Corollary 6.4. For sufficiently large p we have 



Gar = — arccos 

TT 



(l-f)+o(p). 

^ P A 6a — 

Gzi = — arccos 1 + 

^ V 4 y 



as p 00. 
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4. Large 4p-periodic real potentials. Now we consider Shrodinger operator H on 
armchair nanotube with large periodic potentials. We show that in this case the structure 
of the spectrum is the same in the essential as for zigzag nanotube (see Theorem 12.61) . but 
the Proofs are different. 

Theorem 6.5. i) Let v = (fn)lCo ^ 4p-periodic (p > 2) real potential such that Vi ^ vj, 
l^i^ j ^Ap. Let a{t) = a{H\tv)) and auit) = a{Jk{tv)). Then 

N 4p 

fc=l j=l 

where intervals Ckjit) satisfy 

Here Xj are defined in (16.441) and Qi are defined in (16.311) . 

Moreover, if v\ < v\ < .. < and h is sufficiently small, then all intervals (Tk,jif) are 
disjoint for sufficiently large t. 

Proof of Theorem 16.51 Recall that 

a = ak= \ Q l,s = eiv, c/„=( ^^^^^ 1 , n e Z, aa* = h. 

(6.26) 

Also we use notation dn = dn{v), where v = (f i, .., v^p). The monodromy matrix for operator 
Jkiv) is 

Mk = Mk{z) = Mk{z,v)=M2p..Mu Mn=( \ w/' w ^ V (6-27) 

\ O-k a,k\Z an) J 

It is well known that 

a{Jl) = {z : det(Mfc(2) - r) = far some r E §^}. (6.28) 
Using (I6.27P we obtain 

\ f Pi P2 



[ a2p{z - d,p).M^ -d,))"'\p, P4 J ' ^^-^^^ 

where Pj = Pj{z — d2p,--,z — di) is a 2 x 2 matrix polynomial and degPj < 2p for all 
j = 1, ..,4. Also, using (16.261) and periodicity of v, we deduce that 

(6.30) 
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where Qj = Qj{z — v^p, .., z — Vi) are polynomials and degQ^ < 2p, sets Qj are 

Qi = U^;d{4j + 1, 4j + 2}, Q2 = \ Qi = U^Io'{4j + 3, 4j + 4}. (6.31) 
Let Dk{z, r) = Dk{z, r, = det(Mfc - r/2). Using flOOD - flOB we get 

Z}fc(z,T) = r^ + deta^r3( JJ JJ + +det r2(]^(2-t;„) +i?2) (6.32) 

neQi neQ2 n=l 

+r^i + R2, 

where polynomials 

Ri = Ri{z-V4p,..,z-vi), degRi<2p, R2 = R2{z -v^p, z -vi), degi?2 < 4p, (6.33) 

Ri = Ri{z - V4p, ..,z - vi), R2 = R2{z - Vip, ..,z - vi). (6.34) 
are not depended on r. Let t E , z E M., then it is well known, that the polynomial 

D,{z) = D,iz, t) = Dkiz, T, v) = (det a, '^)r-2Dfc(^, t) = 1[{z- v^) + Oiz^P-'), z ^ 00. 

n=l 

(6.35) 

is real, since it has only real zeroes, because the spectrum of Jk is real. Let r G S^, 2; G M, 
then using 06.351) . fl6.32p and Dk{z,T) = Dk{z,T), al = a^^ we deduce that 

^1 = det af{ Y[{z-Vn)+ Y[{z- Vn) +1^1), R2 = det af. (6.36) 
Substituting (16.361) into (I6.32p and using (16.351) we deduce that 

Dk{z) = l[{z-Ur,) + R2 + 2 Re(r det a-^)( Y[ (z - u^) + H {z - t;„)) (6.37) 

n=l nGQi neQ2 

+2 Re(r det a~PRi) + 2 Re{T^ det a'^^). 

Now we denote a = A = | and -Ffc(A) = Ffc(A,a) = Ffc(A,r, a) = t^^PDk{z,T,tv). Then, 
using (16.371) . fl6.33p . we deduce that 

Fk = n(A-i;„) + aGi(A,a)+a2P2Re(rdetafc^)( JJ iz-Vn)+ JJ (^-^^n)) + a'^+'G'2(A, a), 

n.=l neQi neQ2 

(6.38) 

where Gi,G2 are polynomials and Gi is not depended on r. Let Xj{a) = Xj{a,T) be zeroes 
of -Ffc(A) such that Aj(0) = fj, these are analytic functions. Using similar arguments as in 
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"zigzag case", we deduce that derivatives (XjY^^Q) are not depended on r for all j e Nip, 
r G N2P-1 and 



(A,)(2f)(0) = ^2Re(rdeta/)^ ^ ^^^^^ ^ ^ ^^^^ ^ ^ 2. 



These yield 



n 



ne(Qi\j) 



(Vj - Vn) 



(6.39) 



(6.40) 



where j e Qj for some i — 1,2. Let 2;j(t) = Zj{t,T), j e be zeroes of Dk{z,T,tv), then 

4 



Zj — tXj and 



t^p-' n 



+ 0{t-'P), t^oo, (6.41) 



where the spectrum (j{Jk{tv)) — \J^ak,j{t). Introduce the C^p^^p matrices ivfe(r) = Lk{T,t) 
and Bfc(r) = Bk{T,t) by 



Bk + diag(tt') 



where 



f d al ... f\ 

ak d a^, ... 

Ojfc d ... 

y ra^ ... d J 

e^^'s 



+ diag(tt'), 



(6.42) 



(6.43) 



Let Aj(t) = Aj(t,r) be eigenvalues of L^, it is well known, that o'kj{t) = Xj{t,E>^). Then 
perturbation theory gives us 



neN4p\j 

V^-,--l+yj,l+yj,3 _ Sk {Vj,-2Vj,-l Vj,l+Vj, 1 ^^j,2 ^^j,3 ) 



(6.44) 



Vjt 

yl^ _ ^"j,-3+Vj,-l+Vj,l _ Sk{Vi,-3Vi,-2Vj,-l+Vj,-lVj,lVj,2) 



o{t- 



j e 2N 



-t ^ + 0{t--') JG2N + 1 

where Vj,k = (vj+k - Vj)-\ = 2 Re(s'e*(''i+''^-2''^)). 



, i — > 00, 
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7 Appendix 



Below we consider the unperturbed Jacobi operator J°(a) = J(a, 0) given by (see (13.71) ) 
(J°(a)?/)„ = an-iVn-i + dnVn+i, = ct > 0, a2„+i = 1, u e Z, y = {yn)n&- (7.1) 

Lemma 7.1. The eigenvalues z'^ and the eigenvectors e,^, (n, s) G Np x N2 of the matrix 
K^{e^'f', a), G M ( given by (13.101) ) have the forms: 
%fen = a + e-" 7^ 0, r„ = then 

< = (-l)1e.|, e: = (2p)-^(e:,)SiGC2^e;2, = (-l)V^'''", e„,2,+i = e^^-^. (7.2) 

-(f^n = 0, then the eigenvalue = z"^ = has the multiplicity two and the corresponding 
orthogonal eigenvectors are given by 

el = (1, 1, -1, -1, 1, 1, ..)T, el = (1, -1, -1, 1, 1, -1, ..y G C^f. 

Proof. We need the simple fact. Let i^°(r)e = ze for some z,t and the eigenvector 
e = {fnfi ■ Introduce two numbers /o = T~^f2p, /2p+i = t/i- Then 

M2(Z)(/„_1, fnV = (/n+1, /n+2)^, M.izYifo, /l)^ = (/2p, /2p+l)^ = r(/o, /l)^, 

and (/o, /i)^ is the eigenvector of the monodromy matrix M2 given by (14. 2 p at v = 0. 

Conversely, let M2{zi){fo, fi)~^ = r^fo, fi)~^ for some t,zi. We introduce the vectors 

(/n+l, /n+2)^ = M2{z,){fn-l, fnV, U G ^2^-2- Then 

/^"(r, a)ei = ziCi, where Ci = (/„)i^. (7.3) 

Recall that (see (14.31) ) the Lyapunov function A2 (corresponding to M2) is given by A = 
1 Tr 71/2(2;) = ^(z^ — — 1). Using these arguments we will determine the eigenvalues and 
the eigenvectors of the matrix K^{T^a). Firstly, let Zg = Zs{r) be solutions of the equation 
A{z) = cosr for fixed r G M. Then {zg)"^ = a^ + 2acosr + l = e = a + e"', which yields 

^1 = — I £ I , ^2 = I ^ I • 

We will determine the eigenvectors of the monodromy matrix M2{zs), s 
value r = e^^, since A{zs) = cosr. Firstly, if e 7^ 0, then we obtain 

\ —Zs a + 2cosr — e J \ (— 1)^ |£:| 
and the corresponding eigenvectors are given by 



1, 2 for the eigen- 



-l)"kl 



(7.4) 
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Define the vectors = (e^)^^ by 

) = V'{r), ( ^ = MHz,) (I ) = e'^'v'ir). (7.5) 

Then using (17.31) we deduce that K'^{e^P^,a)e'^ = Zgc'^, where identities (17.41) . (17. 5p give the 
components of by 

e^ = (e-)f=i, e^,. = (-l)V^ "^.-i = e'^'''^, 

which yields (17.21) . since solutions of the equation e*^^ = e**^ has the form = ^ + 

Secondly, if e = a + e*^ = 0, then we deduce that a = 1, e**" = —1, Zg = 0, s = 1,2 and the 
matrix M2(-2s) — e^^l2 = 0. The corresponding eigenvectors have the forms r]^ = (—1, 1)^, 
r/^ = (1, l)"*^ and using arguments as above, we obtain the proof of the case e = 0. ■ 

Corollary 7.2. The spectrum of the operator J°(a) given by (17.11) has the form 

a{J%a)) = ull,al cr° = [A+_i, A;], = X^ia) = z^{a,0), X^^ = -X+ = a + l, 

X^ = iy^\a + e'^\, z/^ = (il)''"-'' sign(n - p), n e Ngp-i, sign(O) = 1, 

where Xf^ (and Xf^_^_^) are all eigenvalue of the matrix K^(l,a) (and K^{—l,a)) given by 
(I3.10p . Corresponding eigenvectors o/-ft'°(l,a) (and K^{—l,a)) are given by 

7±^Z^(n\- \ (f^S^^ - iAt^^ _ ±j ±ia.rg(a+T„) 

— ~ X 1 \Jj,n)i=li J2j,n ~ '^n ' n ' J2j+l,n ~ 'n ^ ' 

(2p)2 



Tn = e p , j eNp, a + Tn^O, (7.6) 
and 

z+ = (2p)-i(i, 1, -1, -1, 1, 1, .y, z~ = {2py-2{i, -1, -1, 1, 1, -1, .y, a + x„ = O, 

and A~(a) = A^(a), n G N2p-i \ {p} has multiplicity two. Also A^ (a) < A^ (a), a 1 and 
A~(l) = Ap (1). The vectors and Z^,n G f^2p-i (if^G orthogonal. 

Proof follows from Lemma [TTTl In particular, we have the following identity {2p){Z^, Z~) = 

3D coordinates of r^^ and bj in the case of armchair nanotube. We rewrite similar 
formulas from [BKj adapted for our case 

rn,j,k = {Rcosan,j,k,Rsman,j,k,nh), n G Z, j G {0, 1}, k G Zn, (7.7) 

where 

_ 2n{k-n) _ 27i{k-n) 

C^2n,j,k = 1" CtOj5 0!2n+l,j,k = 1" Ctlj, 
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«o,o = 2p, ao,i = — , ai,o = p - a, = — , 



smQ; = — , sinp = — , ti 



2R' ' K sill;! 



h=^2 + R,R2 - 2i?2, R~=^{jRy-l, j = l,2, 
and the magnetic constants are 

-B(-R2 - -Ri) , BR2 . 
bi=b2 = , 63 = (7. 



References 

[ARZ] Avron, J. E.; Raveh, A.; Zur, B. Adiabatic quantum transport in multiply connected systems. Rev. 

Modern Phys. 60 (1988), no. 4, 873-915 
[BK] Badanin A.; Korotyaev E. Magnetic Schrodinger operators on armchair nanotubes. 

http: / /a rxiv.org/abs/0804.0171 
[Ex] Exner, P. A duality between Schrodinger operators on graphs and certain Jacobi matrices. Ann. 

Inst. H. Poincare Phys. Theor. 66 (1997), no. 4, 359-371. 
[Ha] Harris P. Carbon Nanotubes and Related Structures, Cambridge Univ. Press., Cambridge, 1999. 
[li] lijima S. Helical microtubules of graphitic carbon. Nature. 354(1991), 56-58. 

[Kl] Korotyaev, E. Effective masses for zigzag nanotubes in magnetic fields, Lett. Math. Phys., 83 (2008), 
No 1, 83-95. 

[KKr] Korotyaev, E.; Krasovsky, I. Spectral estimates for periodic Jacobi matrices, Commun. Math. Phys. 
234(2003), 517-532. 

[KKul] Korotyaev, E.; Kutsenko, A. Inverse problem for the discrete ID Schrodinger operator with small 
periodic potentials, Commun. Math. Phys. 261(2006), 673-692. 

[KKu2] Korotyaev, E.; Kutsenko, A. Lyapunov functions for periodic matrix-valued Jacobi operators, will 
be published in " American Mathematical Society Translations" . 

[KKu3] Korotyaev, E.; Kutsenko, A. Armchair nanotubes in magnetic field, tight-binding models, in prepa- 
ration. 

[KL] Korotyaev, E.; Lobanov, I. Schrodinger operators on zigzag periodic graphs. Ann. Henri Poincare 8 
(2007), 1151-1176. 

[KLl] Korotyaev, E.; Lobanov, I. Zigzag periodic nanotube in magnetic field, preprint 2006. 
[KuP] Kuchment,P.; Post, O. On the spectra of carbon nano-structures Commun. Math. Phys. 275(2007), 
805-826. 

[La] Last, Y. On the measure of gaps and spectra for discrete ID Schrodinger operators. Comm. Math. 

Phys. 149 (1992), no. 2, 347-360. 
[vM] van Moerbeke, P. The spectrum of Jacobi matrices. Invent. Math. 37 (1976), no. 1, 45-81. 
[N] D.S. Novikov, Electron properties of carbon nanotubes in a periodic potential. Physical Rev. B 

72(2005), 235428-1-22. 

[Pa] L. Pauling, The diamagnetic anisotropy of aromatic molecules, Journal of Chemical Physics, 4 
(1936), 673-677. 

[Pk] Pankrashkin, K. Spectra of Schrodinger operators on equilateral quantum graphs, Lett. Math. Phys. 
77 (2006) 139-154. 

[RS] M. Reed ; B. Simon. Methods of modern mathematical physics. IV. Analysis of operators. Academic 

Press, New York-London, 1978. 
[RSI] K. Ruedenberg and C. W. Scherr, Free-electron network model for conjugated systems. I. Theory, 

The Journal of Chemical Physics, 21 (1953), 1565-1581. 



28 



[SDD] Saito, R.; Dresselhaus, G.; Dresselhaus, M. Physical properties of carbon nanotubes, Imperial College 
Press, 1998. 

[SI] Simon. B. Orthogonal polynomials on the unit circle. Part 1. Classical theory. American Mathemat- 
ical Society Colloquium Publications, 54, Part 1. American Mathematical Society, Providence, RI, 
2005. 

[S2] Simon, B. Orthogonal polynomials on the unit circle. Part 2. Spectral theory. American Mathemat- 
ical Society Colloquium Publications, 54, Part 2. American Mathematical Society, Providence, RI, 
2005. 

[vMou] van Mouchc, P. Spectral asymptotics of periodic discrete Schrodinger operators. I. Asymptotic Anal. 
11 (1995), no. 3, 263-287. 



29 



